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I. INTRODUCTION 



^ ' Radiative decays cj) rjir^j and 4> ~* ir'^Tv'^'y are studied assuming that these decays are caused 

, through the intermediate ao (980)7 s-nd /o (980)7 states, respectively. Fitting the experimental data 

. of the rjn'^ and tt^tt" invariant mass spectrum in the decays 4> ~* rjir^'y and 7r''7r°7, it is shown 

' that the processes (j> — j- ao7 and (j> ~^ /o7 are dominated by the K^K~ loop interaction rather 

I than the contact </>ao(/o)7 one both for the non-derivative and derivative SPP coupling. The 

r-~r , experimental data of r[(f> foy]/T[(l> aoj] predicts that gfyKR/daoKK ~ 2. Considering the 

' effects of the mixing between low mass scalar qqqq states and high mass scalar qq states to these 

' coupling constants QfgKK and QaoKK^ one suggests that the mixing is rather large. 

(N ■ 
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. For a long time, the radiative decays of the (f) to ir'^rjj and 7r°7r°7 have been analyzed, assuming that the decays 
■ (/) — > 7r°77(7r°}7_proceed through the (j) — > ao(/o)7 decays, to reveal the structure of the light scalar mesons ao(980) and 
. /o(980) [il, These analyses are performed assuming the charged K^K~ loop diagram in the coupling (jjo^oifo)"/ 
O i' and the contact (vector dominance) (/)ao(/o)7 interaction. In this work, we analyze the data for the ir^r] and tt^tt^ 
I invariant mass spectrum in — > 77^777 and — > 7r''7r''7 decays given in recent precise experiment Q, assuming charged 
K'^K~ loop diagram and contact 0/o(ao)7 interaction. For the S'(scalar meson)-P(pseudo-scalar meson)-P(pseudo- 
(-H I scalar meson) interactions appeared in these decaying processes, we consider the cases non-derivative interaction and 
• • • derivative one, latter of which is adopted in the literature [3|. The result obtained in our present analysis shows that 
these processes are caused through the K~ loop diagram dominantly. 

When the K^K^ loop diagram is dominant in the decays (j) — > /07 tt^tt^^ and 4> — > ao7 7r"r/7, the value of the 
ratio g fgKK / 9aoKK is obtained from the experimental ratio T{<p aoj)/T{(j) — > /07), that is the rather large value 
9 foKK 1 9aoKK ~ 2. These coupling constant strengths depend on the structure of the scalar mesons, that is, these 
scalar mesons are constituted of qq or qqqq, or are mixing states of qq and qqqq. Authors in Ref. ^ argue that the 
data of ^ o-o{fo)"l ■K^r]{-n^^ decays gives evidence in favor of the qqqq nature for the scalar ao(980) and /o(980) 
mesons, and authors in Ref. argue the matter of mixing between qq and qqqq states. 

Recent many literature (refer to the "Note on scalar mesons" in fs^) suggest that the low mass scalar nonet (/o(600), 
is:o*(800), ao(980), /o(980)) are qqqq state and the high mass scalar mesons (ao(1450), i^o*(1430), /o(1370), /o(1500), 
/o(1710)) are conventional L = 1 qq nonet plus one glueball. We assume a strong mixing between low mass and 
high mass scalar mesons to explain the fact that the high L — \ qq scalar nonet are so high compared to other 
L = I qq l"*""*" and 2++ mesons ^iTj. Assuming that the coupling strengths causing the mixing between I — 1 ao(980) 
and ao(1450), / = 1/2 K^{800) and i^^(1430) and / = (/o(600), /o(980)) and (/o(1370), /o(1500), /o(1710)) are 
same, we analyzed the S — > PP decays using derivative SPP couplings <§\. Fitting the various experimental SPP 
decay widths, we obtained the mixing angle between ao(980) and ao(1450) as ^^9°. In our previous work Q, we 
analyzed the T((j) —^ /07) and r(0 007) assuming the contact (vector dominance) coupling for the ao(t>j and fo4>"t 
interaction and using the mixing strength obtained in previous work and then suggested the importance of the 
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mixing effect for the explanation of the rather large ratio T((j) — > /o7)/r(0 ao7)- 

In section 2, we analyze the data for Tr^r? and 7r°7r° invariant mass spectrum of the dBR{(j> ir'^ ir'^ ^) / dq and 
dBR{(j> ■K^'q^)/dq assuming the intermediate scalar states /o(980) and ao(980). In this analyses, we consider the 
contact and K~^K~ loop interaction for 0/0(00)7 coupling, in cases of derivative and non-derivative SPP coupling. 
In section 3, we reanalyze our mass formula for low mass nonet scalar and high mass nonet scalar + glueball adopting 
the new mass data of ifpl^OO) [1]. In section 4, we express the SPP couphng constants gaoir-K, S/otttt etc. using the 
mixing parameters between low and high mass scalar mesons. We pursue the best fit analyses for the S — > PP decay 
data using the mixing parameters obtained in section 3 for both non-derivative and derivative SPP interactions, and 
then obtain the best fit gf^^R etc. Comparing the best fit gj^KK etc. with the values obtained from the <j} — > fo{aQ)^ 
decays, we suggest that the non-derivative coupling is reasonable than the derivative one and the mixing between qq 
and qqqq is rather large. 

II. ANALYSIS OF THE <j) n°r]-f AND ^ 7rV°7 DECAYS 

A. — * TT^ri-y decay 

Firstly, we consider the decay (j) —^ ao(980)7 — > ir^-qj shown in Fig. 1. The invariant mass distribution of the 
branching ratio dBR{(f> ao (980)7 ^ ■n^'q^)/dm is expressed as (refer to first paper in Ref. [l|) 



dBR{(j) — !■ ao7 n'^rjj) 2vr? 1 r(0 — s- ao7 : ra)V{aQ Tr'^r/ : m) 



(1) 



dm TT \Dao{m'^)\'^ 

where F^ is a decay width of <p and l/I?Q„(m^) represents the propagator of the intermediate state oo, 

Dagim^) = rr? - ml - inia^TaQ. (2) 



F(ao — > TT^rj : m) is the decay width on the virtual mass m of intermediate ao defined as m — \/ Qo <i^, 



F(ao TT rj : m) 



Snm'^ 2m 

1 for non-derivative coupling, 



2 



m^-ml~ m2 \ , , . ^. (3) 
lor derivative coupling, 



where coupling constant gaam] is defined as 



7i,r/ / N 0/ \ . / w I 1 for non-derivative coupling, ,.s 

M(ao(g) ^ 7: (qi) + r^iq,)) = \ q, . q, for derivative coupling. 



K 



FIG. 2: Diagrams for the contact and K'^K loop coupling contributing to g4,aoi{'m) 
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FIG. 3: Diagrams for the K loop coupling contributing to gtpaojin^) 



T{(j) ^ ao7 : m) is the decay width on the virtual mass m — -y/^Q — of intermediate state qq, 

/ 2 ? \ 

Q, / rrj, — m \ 



2m^ J 

where coupling constant grtiao-iifn) is defined as 



(5) 



(6) 



For the coupling (70QQ^(m), contact interaction and K^K loop interaction contribute as shown in Fig. 2, and then 
gtpaa-yi^) expressed as 



(7) 



9^i^j°°^{Tn) is calculated for non-derivative gqK^K coupling by many authors (N. N. Achasov et al. and F. E. 
Close et al. in [ij) considering three diagrams (a), (b) and (c) shown in Fig. 3, as 



KKloor,, X 94,KK9aoKK rl u\ c j ■ r 

3^007 ~ 2 ■ — 2 — ^(ajb). tor non-derivative couplmg 



(8) 
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The quantities a, b are defined as a = m'^/rrv^, b = m? /m?j^ and /(a. 6) arisen from the loop integral is 

1 2 



/(a, 6) 



where 



2(a - 6) (a - 5)2 



/ 



/ 



2V^ 



(a - 6)2 



a; > 



(9) 



log ^ ITT 



X < 



\/Ax - 1 sin ^ 
1 



Vl -4a; 



2V^ 



log ZTT 



X > 



X < 



1 

2x 



(1 ± Vl - 4x) 



The coupling constant g^KK is defined as 

M(</.(p, e^) ^ K+{q,) + i^-fe)) - S^i^xe^l?^ ^ 92^, 
and decay width is expressed as 

„2 „ / 2 \ 3/2 



rf 



K+ + K- 



An 3m2 I 4 

The coupling constant g^KK is estimated using the experimental data r((^ — > + if 

g^^^^f =4.55 ±0.06. 

For the a^K^ coupling, ga„KK is defined by the similar expression as Eq. (4) 

n r/ I \ T^+t \ , T^-t \\ f 1 for non-derivative coupling, 

M{a,{q) ^ (gi) + X (g^)) - ff.^KK >< | . for derivative coupling. 



(10) 



(11) 



(12) 

2.10±0.05MeV Q as 

(13) 



(4') 



For derivative coupling of a^K^K , K^K loop diagram contribution 9^^^^°°^ {fn) is calculated by D. Black et al. 
[2| considering four diagrams (a), (b), (c) and (d) shown in Fig. 3, as 



KK loop / \ 



9<j,KK9aoKK 2mj^ - 
27r2im|. 2 



/(a, 6) for derivative coupling. 



(8') 



Using Eqs. (2), (3), (5), (7), (8), (8'), we parameterize Eq. (1) as 
dBR{(j) ao7 71^777) 



dm 



= Gi- 



|G2 + i 


2m^ — 
2 


/(a,6)|2 


\G2 + \ 


2m|^ — rn^ 


/(a,6o)P 


2 



9 2 

— m 



m, 
m 



a 



/ (m2 — (m,, + m7r)2)(TO2 — (to^ — m7r)2) 
(m2 — m2)2 -|- m2r2 V (m2 — (m,, + m.„)'^){m'^ — {171^^ — niTr)^) ' 



(14) 



where Gi, 62, 60 ^^'^ defined as 



Gi - 



rFf 



7rr0r2 

^contact / / 94iKK9ao 

^2 — g^fa / 2 



ao7 : ma)r(ao ^ 7777° : m^), 



if 



(15) 



60 



TO' 



TO 
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FIG. 4: dBR{4> -k*^ r)^) / dm (in unit of GeV~^) as a function of the -k^-t) invariant mass m (in GeV~^). Solid line shows 
the best fitted curve for the non-derivative SPP couphng interaction and the dashed line shows the best fitted curve for the 
derivative one. Experimental data indicated by circles are from the SND collaboration in Ref. and those by filled circles 
and filled squares are from KLEO collaboration in Ref. [J]. 



are replaced to 1 for non-derivative SPP coupling. r(ao — *■ rjir^^ma) and 

T{(j) —f aoj,ma) are defined in Eq. (3) and (5) settling m ma- We fit the Eq. (14) varying the parameters Gi and 
G2 using the experimental data from the SND collaboration and KLEO collaboration in Ref. [J . Best-fitted curves are 
shown in Fig. 4; solid line for non-derivative SPP coupling and dashed line for derivative SPP coupling are obtained 
for the choice of the parameters Gi and G2 as Gi — 4.1 x 10~^GeV~^, G2 = —0.16 for non-derivative coupling and 
Gi — 3.9 X lO^^GeV^^, G2 = 0.08 for derivative coupHng. For these choices, the estimated BR{(j) — > ir^r]^) are 
estimated as 

Gi = 4.1 X 10-4GcV"\ G2 = -0.16, BR{(l3 tt^t]-/) = 7.03 x 10"^ 

for non-derivative coupling, (16) 
Gi = 3.9 X 10-4GeV~\ G2 ^ 0.08, Bi?((/i ^ 7r°777) = 7.12 x 10"^ 

for derivative coupling. (16') 



and factors 



2ml 



and 



Estimated value for BR{(j) — > ir^rj^) is consistent with the experimental data BR^^P{(j> — > 7r°?77) = (8.3 ±0.5) x 10^^ 
[3, Q. Also, estimated value of Gi is consistent with the value evaluated from the experimental data (Ref. [5]) using 
the relation Eq. (15), 

Gf P = -^BW'P{(j) ao-f : ma)BR''''P{ao t]t:° : nia) = (5.96 ± 2.47) x lO-^GcV'^ 

TtI a 

Furthermore G2 is very small compared to \I{a,bo)\ = 0.902 for iria = 0.985 GeV, then we can suppose that the 

loop contribution is dominant in the (j) ~^ t^^t]^ decay. Supposing that the decay (j) ~^ ao7 is caused through 
only the K~ loop interaction, we obtain the result. 



ao7) = 



9^KK9aoKK 



2m\ 



/(a, 60) 



2m0 



(17) 



where the factor 



2m\ 



is replaced to 1 for the non-derivative coupling. Using the value Eq. (13) of g^KK 
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and the experimental value r{(j) ao7) = (0.323 ± 0.029) x lO'^MeV in Ref.Q, we obtain the result 



_ 2.18 ±0.12 GeV, for non-derivative coupling 
ffaoKi? = i 9.04 ±0.50 GeV-^ for derivative coupling 



(18) 



Using relations Eqs. (15), (3), (17) and estimated results Eq. (16), (16'), (18), we obtained the values for gaoTrrj, 

_ J 1.89 ±0.75 GeV, for non-derivative coupling 



B. 



7r''7r'^7 decay 



For the decay (j) foj tt'^tt'^j, the invariant mass distribution of the branching ratio dBR{(j> — > foj — > Ti^TT^^)/dr, 
is expressed similar to Eq. (14) for the case — > ao7 — > ■k'^ti^ as 

dBR{(j}^ fo-/ TT^TT^) 



dm 



= Gi- 



\G2 + \ 


2m^ — 
2 


/(a,6)|2 


\G2 + } 


2 


/(a,6o)P 


/ 


— 4mJ 





TO 



(m^ — TOy)2 ^ TO^Fj y TO^ — 4toJ ' 



where Gi, G2, 6o are defined as 



Gi 



7rF^F2 



/o7 : m/)F(/o ^ tt tt : m/), 



"-^2 — 94, 



contact / / 94>KK9foKK 
'^''■f ' \ 27r2TO2 



K 



bo 



Here, F(/o 7r°7r° : m/) and F(0 /07 : to/) are expressed as 



r(/o ^ : m/) 



TOj — 4to2 



167rm2. 



TO^ — 2m2 



for non-derivative coupling, 

2 

for derivative coupling. 



where coupling constant g/o,r7r is defined as 



M(/o(<z) ^ ^°(gi) + ^°(<72)) = 7r.g/o.. X 



1 for non-derivative coupling, 
q\ ■ q2 for derivative coupling, 



and 



r((?!)-> /o7 : »7i/) 



.2 _^ I ~"^/ 



r50/o7(™/) 



2TOrf 



2to^ — rrij 



/(a, 60) 



In Eq. (24), factor 
as Eq. (22), 



2to|- — TOj 



(19) 



(20) 



(21) 

(22) 

(23) 
(24) 



is replaced to 1 for non-derivative coupling. 9f„KK is defined in the similar equation 

(22') 



nr^4^/ \ T^+i \ , T^-i \\ J 1 for non-derivative coupling, 

M(/o(9) - K (gi) + K (q,)) = g^^^^ ^ \ . for derivative coupling. 
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FIG. 5: dBR{<f) ■n'^-K'^'y) / dm (in unit of GeV~^) as a function of the tt^-tt'^ invariant mass m (in GeV~^). Solid line shows 
the best fitted curve for the non-derivative coupling and the dashed line shows the best fitted curve for the derivative one. 
Experimental data indicated by circles are from the SND collaboration in Ref. [j, and those by filled circles are from KLEO 
collaboration in Ref. Q. 



We fit the Eq. (19) using the experimental data from the SND and KLEO collaboration in Ref. and obtain 
the best-fitted curves shown in Fig. 5; solid line for non-derivative coupling and dashed line for derivative one. The 
choice of the parameters Gi and G2 for these best fit and estimated value for BR{(f) — > n'^rij) are obtained as 



Gi = 7.1 X 10-''GeV"\ G2 = 0.001, BR{(j) -> 7r"7r"7) = 1.06 x 10""* 

for non-derivative coupling, 

Gi = 6.9 X 10-'*GeV"\ G2 = 0.055, BR{(j3 ttV^) = 1.08 x 10""* 

for derivative coupling. 



(25) 
(25') 



Estimated value for BR^cj) — > 7r°7r°7) is consistent with the experimental data BR°^'P{(j) ~^ 7r'^7r°7) — (1.09±0.06) x 10 
[3,01, and estimated value of Gi is consistent with the value evaluated from the experimental data (Ref. [^), 



Gf P = _Bi?«'^P( 
ttE/ 



/o7 : mf)BR-^^(f^ 



TT TT 



TO/) = (10.0 ±4.8) X lO^^'GeV" 



As the case for the decay t^'^VIj contact 50/07 interaction (G2) is very small compared to |/(a,5o)| = 0.783 for 
ruf — 0.980 GeV, then one can suppose that the K^K^ loop contribution is dominant in the </> tt^tt^^ decay. 
We suppose the (j) /07 decay is caused from the K^K^ loop interaction, then can estimate the coupling constant 
9aoKK from the relation 



r(0 ^ /07) 



9(l>KK9foKK 



2m\ 



I{a,bo) 



2md 



(26) 



where the factor 



2to^ 



is replaced to 1 for the non-derivative coupling. Using the value Eq. (13) of g^KK 
^ fg-f) = (0.323 ± 0.029) X lO^^^gy in Ref.Q, we obtain the result 

iTr\l m rr 

(27) 

Using relations Eqs. (20), (21), (26) and estimated results (25), (25'), (27), we obtained the values for Qf^n-n^ 



and the experimental value r(0 — > 

9foKK 



4.72 ± 0.82 GeV, for non-derivative coupling 
20.0 ±0.50 GeV^^. for derivative coupling 



1.12 ±0.69 GeV, for non-derivative coupling 
2.43 ±1.50 GeV~^. for derivative coupling 



(27') 
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Fig. 6 OZI rule allowed graph for qqqq and qq states transition 



The rather large value of the ratio g f^KK / 9aoKK ^ 2 suggests that the ao and fo scalar mesons are not the pure 
qqqq states but there exist the mixing (inter-mixing) between qqqq and qq scalar mesons. Furthermore, the existence 
of the coupling g/oTru- suggest the intra mixing between qqqq /o(600) and /o(980) scalar mesons. 



III. MIXING BETWEEN LOW AND HIGH MASS SCALAR MESONS 



In this section, we review the mixing among the low mass scalar, high mass scalar and glueball discussed in our 
previous work 0, The qqqq scalar SU{3) nonet S'^ are represented by the quark triplet qa and anti-quark triplet 
q"^ as 

~ e^'^q^qaebeffq^ (28) 
and have the following flavor conflguration [13], 0: 

dssu, -^{dsds — susu) , suds <;==^ aj, Aq, 

dsud, suud, udsu, udds -4=> , kP , n'^ , k~ 

-^{dsds + susu) ^=4> Jns ~ /o(980) 

udud /nn ~ /o(600) 

The high mass scalar mesons 5'^° are the ordinary SU (3) nonet 

Ob ^ q qb- 

The inter-mixing between qqqq and qq states may be large, because the transition between qqqq and qq states is 
caused by the OZI rule allowed diagram shown in fig. 6. Considering the above flavor configuration for qqqq states, 
the expression for this transition is suggested as 

Lint - Aoi [a+a'- + a-a'+ + aj]a'° + K*+K'*- + K*' K'*+ + A'*°K^*° + KfR'f 

+V2fNNfN + Ins In + ^fNsf's\- (29) 

The parameter Aqi represents the strength of t he inter -mixi ng and c an be considered as rather large. When we 
represent the I ~ \ pure qqqq and qq states by ao(980) and ao(1450), and masses for these states by m ? ^^^^^ and 

m?-^ -, the mass matrix is represented as 

ao(1450) 



ao(980) 
•^01 



• (30) 

'00(1450) / 

Diagonalising this mass matrix, we can get the masses for the physical states ao(980) and ao(1450) represented as 
mixing states of ao(980) and ao(1450); 



ao(980) = cos g, ao(980) ~ sin 0a ao(145O), . 
ao(1450) = sin6laao(980) -f cos6laao(1450). ^ ' 

Mixing angle 6a and before mixing state masses ttz^^ (980) (i450) represented by the inter-mixing parameter 

Aoi as 



2 2 

^ao(1450) ~ ^ao(980) 



'2 2 \ *^ 

"^ao(1450) ~ "^ao(980) j _ ^2 



3 a = tan-1 (32) 



ao ( 1450) ~ V "^ao( 1450) + ^a, "^ao (980) ~ V "^1,(980) 
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Fig. 7 OZI rule suppression graph for qq — qq transition. 




where ^Ti^^i-ggo) and TOt„g(i45o) are the masses of the states ao(980) and ao(1450). 

Similarly, for the / = 1/2 Kq{800) and Kq{1A30) mesons, the mass matrix is represented as 

/ ^2 ^ 



K*(800) 

Aoi 



01 



(33) 



if* (1430) 



where 



'ifjCsoo) 



and 



'if* (1430) 



are the masses of pure qqqq and qq states Kq{800) and Kq{143()). The physical states 



^^0(800) and Kq{1430) are written by the before mixing states Kq{800) and ii'o(1430) and mixing angle 6k as 



if* (800) = cos6>k Xo*(800) - sin (1430), 
ii'*(1430) = sin6»Ki^o(800) + cos eix-fi'o (1430). 



(34) 



Mixing angle 9k and before mixing state masses itI k' (800) ^^'-^ ^ k'{i430) represented by the inter-mixing parameter 
Aoi as 



( 1430) (800) 

= 2 



Ok = tan - — 

Aoi 



i 



™/s:*(i430) ™/s:*(800) 



- >2 
'^Ol' 



(35) 



"^i£''(1430) ~ ■\/"^^*(1430)"+^' '^A:J(800) ~ ■\/^"^*(800) 



where tokj(800) and ^^^.(1430) are the masses of the physical states i4rg(800) and Kq{IA3Q). 

Next, we consider the mixing between 7 = low and high mass scalar mesons. Among the / = 0, L = 1 qq scalar 
mesons, there are the intra-mixing weaker than the inter-mixing, caused from the transition between themselves 
represented by the OZI rule suppression graph shown in Fig. 7, and furthermore the mixing between the qq scalar 
meson and the ghieball caused from the transition represented by the graph shown in Fig. 8(a). Thus, the mass 
matrix for these I = 0, L = 1 qq scalar mesons and glueball is represented as 



m%, + 2Ai 
y2Ai 
V2Xg 



m 



N' 



m- 




(36) 
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Fig. 9. OZI suppression graph for qqqq — qqqq transition. 



where Ai is the transition strength among the I — qq mesons, A^ is the transition strength between qq and gluebaU 
gg and Xgg is the pure glueball mass square. For the hght / = qqqq scalar mesons, there are the intra-mixing 
caused from the transition between themselves represented by the OZI rule suppression graph shown in Fig. 9, and 
the mass matrix for these / = qqqq scalar meson is represented as 



+ Ao 



V2Xo 



V2Xo 



(37) 



TO 



NN 



2m4 



TO 



NS 



ao(980)' 



K^(800) ao(980)' 

where Aq represents the transition strength between / = qqqq mesons. 

The inter- and intra-mixing among / = low mass and high mass scalar mesons and glueball is expressed by the 
overall mixing mass matrix as 



(38) 



Diagonalising this mass matrix, we obtain the eigenvalues of low mass and high mass scalar mesons 7 = states 
/o(600), /o(980), /o(1370), /o(1500) and /o(1710). The eigenstates of these scalar mesons are represented as follows; 



/ 
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\/2Aoi 
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V2Xo 


m%g + 2Ao 


Aoi 


\/2Aoi 
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I Rfo{600)NN 



[Rfo{M)l] - 



R 



R 
R 
\R 



'foi980)NN 



( Inn \ 

Ins 
In' 
fs' 
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Rfo{600)NS 
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(39) 



R 



'/o(600)JV' ^/o(600)S' ^/o(600)G \ 



/o(1370)ArAr -n-/o(1370)ArS Rfo{1370)N' -n-/o(1370)S' -n-/o(1370)G 



■/o(980)JV' Rfo{980)S' R, 
fnCi.iynw Rfr,(^R7^l^\s> Rf 



■/o(980)G 



/o(1500)AfAr ^/o(1500)AfS -n-/o(1500)JV' ^fo{1500)S' ^fo(1500)G 



R 



Rf 



Rf 



fo{1710)NN Rfo{1710)NS -"-/o (1710)7V' -n-/o(1710)S' -n-/o(1710)G 



R 



Rf 



Rf 



1 



Using the inter-mixing parameter Aqi and the mass values, TO^^^jg^o) = (0.9848 ± 0.0012)GeV, maj,(i45o) — (1.474 ± 
0.019)GeV, TO^.(8oo) = (0.841 ± 0.030)GeV and to^j(i43o) = (1.414 ± 0.006)GeV, we obtained the mixing angles 6'a, 



9k and before mixing states masses m 



ao(980)' ™ao(1450)' "''if5(800)' '"K* (1430) 



TOi 



TOi 



from the relations (32) and (35). Using 



TO 



vnrj^j^ and rrq^g obtained from the second equations in Eqs. (36) and (37), and parameters 



the values of to^/, 

Ao, Ai, Ag and AgG: we diagonalize the mass matrix Eq(38). When we fit the eigenvalues obtained to the following 
experimental mass values [5|, 



w/o(60o) = 0.80 ± 0.40GeV, TO/„(g8o) = 0.980 ± O.OlOGeV, 

to/„(i37o) = 1.350 ± 0.150GeV, to/„(i5oo) = 1.507 ± 0.005GeV, (40) 
m/o(i7io) = 1.718 ±0.006GeV, 

we obtain the allowed values for Ag, Ai, Ag and Agg- We tabulated the Qa-, 9k, Ag, Ai, Ag, Agg? and Rfg(Q80)NN, 



R 



fo{980)NS, Rfo{980)N', ^fo{980)S', ^fo{980)G 



Rf 



Rf 



for the various values of Agi in the Table I. 



IV. COUPLING CONSTANT gspp AND MIXING BETWEEN qqqq AND qq SCALAR MESONS 



In this section, we first express the (?spp's by the mixing angle 9a, 9k and mixing parameters R/oNS 6tc. Next, 
we obtain the values of the gspp using the various S — > PP decay widths and compare these values with the ones 
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Aoi(GgV2) 




V ) 






-^0 v^*^ * ) 


^(^Lxy^'^^ } 




-^/o(980)A^S 


-R/q(980)A^' 


%o(9SO)S' 


-^/o(980)G 




0.20 


9.7 di 0.5 


9.0 ± 0.5 


n ni Q _|_ n OOQ 


n 07c; -L n 007 


04 -k 04 


/I 1 cro -L n oo8^2 

y L-LO^i ^ u.uuo^ 




—0 Q79 -k n 009 
— yj.iyl^ ^ u.uuz 


n n«c: _l n OOfi 


n OOfi -L n 004 

U..ZZU ^ U.UUI 


_0 01 -k 01 

— U.U-LU ^ U.U-LU 




0.25 


IZ .O ^ u.u 


11 4 -k 


n 0*^9 -1-0 010 


9fi4 -k 008 


0^1 -k 0^1 
u.Uu zn. u.Uu 


/I c:i -L n 007\2 
^^I.OIZ ^ U.UU( j 


—0 097 n C\9f{ 
— U.UZ I ^ u.uzu 


—0 Q'^A OO'i 


n nof. _|_ n 008 


984 -U OO'^ 

U..ZO'± ^ U.UUtJ 


_0 01 R 01 R 

— U.UXU ^ U.U-LU 




0.30 


1 c: n _|_ n « 


1 Q Q _L n Q 

±o . o u .o 


n ncrn -1- OOQ 
u.uou ^ u.uuy 


9c;9 -U OOQ 
yj.i^'j^ ^ u.uuy 


04 04 


/I CI _i_ n 008 

yi. .0 ^ U.UUO^ 


—0 04fi 094 


—0 004 


n 1 1 n -1- OOQ 


"^41 -U OOfi 
U.O'iX ^ u.uuu 


_0 01 fi 01 R 

— U.UXU ^ U.U-LU 




0.35 


1 7 R -1- 1 


1 4 -k 1 


n 079 -1-0 019 


9Q'i -U 008 
u.zoo ^ u.uuo 


0^1 -k 0^1 
u.uu zn. u.uu 


/I r:i 1 _L n 008^2 
^l.O-L-L ^ U.UUO^ 


\J.\J\JiJ _l_ L/WiitJ 


—0 Q09 + n 007 


f) 1 Af) _|_ n ni 2 

W.-LrtW _1_ U.L/-Lii 


n 401 + 007 


—0 094 + 094 




0.40 


20.8 lb 1.2 


19.1 ± 1.2 


0.104 ±0.012 


0.213 ±0.009 


0.05 ± 0.05 


(1.509 ± 0.006)2 


-0.094 ±0.021 


-0.864 ±0.010 


0.178 ±0.014 


0.461 ±0.007 


-0.028 ±0.028 




0.45 


24.2 lb 1.6 


22.1 ± 1.5 


0.146 ±0.014 


0.178 ±0.007 


0.04 ± 0.04 


(1.506 ± 0.002)2 


-0.116 ±0.021 


-0.813 ±0.011 


0.226 ±0.015 


0.523 ±0.006 


-0.014 ±0.014 





TABLE I: The values of mixing angles da, Ok, and the transition parameters Ao, Ai, Ag Agg, and mixing parameters Rfa{9SQ)NN , 
Rfo(980)NS, Rfo{980)N', Rfo(980)S' , Rfo(980)G for the various values of Aoi. 



S 




Fig 10. SPP, S'PP and GPP coupling 



obtained from (/) decay, and then suggest the importance of the mixing between qqqq and qq scalar mesons. 

We use the following expressions for S{qqqq scalar mesonVPP, S'{qq scalar meson)PP and G'(pure glueball)PP 
coupling with coupling constants A, A' and A" , respectively @, 

Lj = Ae-'^eaesStP;^P!+A'S'^{P^, Pn + A"G{Pl P^}, 

for non-derivative coupling (41) 
Ae^'^'^ea.jStd^P^d^Pi + ^S'^{d^P^, d^P^ + A"G{d^Pl d^P^- 

for derivative coupling (41') 

These interactions are represented graphically by the diagrams shown in Fig. 10. We define the coupling constants 
9 SPP' in the following expression, 

V + 9aa-iTr]' O'0'[d'^]TT[dfi]ri + 9a'gTrri'(J-'o' 

+ 9KSKAmKT-[d^]7TK* + H.C.) + gKi'KA[d'']KT.[d^]7,K*' + H.C.) 

+ gK',Kn(Kl[d''\K[d^r, + H.C.) + gA'*'x,(<[a^]i^ [a^]r/ + H.C.) 

+ gK'.Kn'iK^ld^md^]!^' + H.C.) + gK:,'Kn'(K^[d'']K[d,]Ti' + H.C), 

+ gMM)..\h{M)[d^^]iz\d,,]7T ^ gf^^j^,^^T^UM)[^^^^ 

+ 5/o(M)w/o(M)[5l77[ap]'/ +5/o(M),','/o(M)[air7'[a^],7', (42) 

where [9^]'s are replaced to 1 for non-derivative couplings. These definitions of gspp^s are the same as ones of jspp 
in our previous work [i| except for ga^KK^ 9k*K7t, which are related as V^gagKK = laoKK: V^gx^Kyr = Jk^Ktv- 
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Then the couphng constants gspp's are expressed as 

9ao{980)KK = \/2(^ COS 6»a - A'sin(9„), 

9ao{1450)KK = \/2(^ sin (9^ + A' COS 61a), 

9ao{980)nrt = 2(Acos6ia sin^p - \/2A' sinOa cosOp), 
500(1450)707 = 2{A sin Oa sin Op + \PlA! cos Qa COS 6'p), 
fl'ao(i450)7r77' = 2(- A sin Qg. COS + \/2A' COS sin 6lp), 
fl'if5(800)7rif = \/2(^coseK - A' sin 61a:), 

5if5(i430)7rx = ^(Asin6'K + A'cos6'/f), (43) 

— 2{-ARf^(^M)NN + V^^'-R/o(M)JV' + 2A"i?j(,(M)G), 
9fa(M)KK = V^(-^-R/o(M)JVS + ^'RfoiM)N' + V2^'%o(M)S' + 2V^A"i?/jj(M)G)) 

9fo(M)m = 2 ^-^-R/o(M)Jvs cos^p sin^p + iAi?j^(M)jvjv cos^ Op 

+ R}o(M)N' cos^ 6'p + A'Rf^(M)S' sin^ 6'p + ^"-R/o(m)g^ , 

9fo{M)r]r]' = 2 ^^-R/o(M)Ars COS 2^p + -ARf^(M)NN sin 2^p 

+ -^^'-R/o(M)JV' sin26ip - ^'%(M)S' sin26ip^ , 

9fo{M)n'ri' = 2 ^A%(M)ivscos6'psin6'p + ^A%(m)jvjv sin^ 6'p 

+ -^^'-R/o(M)Ar' Sin^ 0P + A'RfoiM)S' COS^ 6'p + A"Rfg(^M)G^ , 

where 6p is the jj-ry' mixing angle related to the traditional octet-singlet mixing angle ^o-8 as 6p = 9o-s + 54.7°. 
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Decay widths for these scalar mesons are expressed by using the couphng constant gspp as 



r(ao(M) ^ K+{mi)K-{m2) + i^°(mi)X0(m2)) = 2 



9ao{M)KK \q\ 

ao(Af) 



' Mm iT?x 2 



r(ao(M) ^7r(TOi) +77(7712)) 



r(ao(Af) ^^(77li)+77'(77l2)) 



9ao(M)'Kr, \q\ 

Stt 777^ 

2 

9aa{M)7rrj' \q\ 

ao(i\i) 



T{K*+{M) 7r+(771l)if°(77l2) + 7r°(777l)if+(7772)) = 3 



9K*(M)nK |q| 



8^ '^K5(M) L 



2 ^ 2n 

'^Mmim-2 



TifoiM) ^+(7ni)^-(777 2) + 7r"(777l)7r'^(7772)) = - 



3 9fo(M)^ 



2 Stt 777^ ,.,s 



r(/o(Af) ^ ii:+(777l)K-(7n2) + X°(777l)X"(7772)) = 2 



r(/o(Af)-> 77(7771) +7?(7772)) = 2 



r(/o(M)^ 77(7771) + 77' (7772)) 



2 

9fo{M)nn' \q\ 



r(/o(Af)->77'(777l)+77'(7772))=2 



°" "Vo(M) 

9foiM)r,'n' |g| 



A'/mim2 



(44) 



fo(M)KK \q\ 
8^ "^/o(M) 



Mmim2 



2 N 2n 
2 

-,2 

A/mim2 



m 



Mm\m2 



Here and rnMmim2 defined as 



\q\ = 



/ Af 2 + — 777^ 
V 2M 



"7Mmim2 = ^yJVP -mf-rr 
and for the case Af w n7i + 7772, we use the next formula for |q| [Tlj |. 



|q| =Rc 



27rr/, 



Af+oo 



M JM-00 




— m\ dm, 



(45) 



where Tm is the decay width of particle with mass M . 

We use the experimental data of the scalar meson decays cited in Ref. [s^; r(ao(980) 7r77 + KK) = 75 ± 25MeV, 
r(ao(1450) ^ 7r77 + 7r77' + KK) = 265 ± 13MeV, r(if*(1430) ^ ttK) = 270 ± 43MeV llJ],_r(/o(980) ^ tttt + KK) = 
70 ± 30MeV, r(/o(1370) ^775 + KK + 7777) = 214 ± 120MeV (H, r(/o(1500) ^ tttt + i^if + 7777 + 7777') = 55 ± 9MeV 
fli ]. r(/o(1710) TTTT + if/C + 7777) = 137 ± 8MeV, for the best fitting of our model parameters, A, A', A" and 
Op. Results are tabulated in Table II for non-derivative coupling case and in Table III for derivative coupling 
case. For the Op, we search the best fit value in the range (54.7 ± 18)°. We estimate the best fit values of these 
parameters for various points (0.2, 0.25, 0.30, 0.35, 0.40, 0.45)GeV^ of inter-mixing parameter Aqi. We show the values 
of r(ao(980) — !■ 7r77 -I- KK) etc. and ga^9%o)-m) etc. for best fitted A, A' , A" and Op. Values in the row just below the 
one denoting r(ao(980) nrj + KK) etc. denote the experimental values and the values in the row just below the one 
denoting (7ao(980)ir?) etc. are the values Eqs. (18), (18'), (27), (27') obtained from (j> aoj/n'^rij and /qj/tt'^tt'^^ 
decay analyses. 



These results in Table II (non-derivative coupling case) show that the values of g, 



'ao(980)7r77 



etc. obtained for Aqi = 



0.30 ~ 0.35(GeV^) are close to the values obtained in — > aoj/n'^ri'j and Joi/t^^t^^I decay analyses. For 
the derivative coupling case (showed in Table III), the values of gao(980)7rr; etc. except for gfg(gso)KK obtained for 
Aqi — 0.30 ~ 0.35(GeV^) are also close to the values obtained in ^ aQj/n^rj-f and — > /o7/7r°7r°7 decay analyses. 



But characteristic feature g^^ 



o{980)Kk/ 9aoi980)KK 



2 obtained in </> — > aQj/n'^ru and — s- /oj/tt'^tt^J decay analyses 



cannot be taken in any values of Aqi for derivative coupling case. 
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Am (GeV^) 


A 


A' 


A" 


dp 


p 

ao(980)^7rT7 + KK 


GeV 


GcV 


GcV 




0.075 ± 0.025GeV 


0.20 


2.8 


1.2 


—0.25 


36.7 


0.189 


0.25 


2.5 


1.2 


—0.23 


36.7 


0.133 


0.30 


2.3 


1.2 


-0.24 


36.7 


0.097 


0.35 


1.9 


1.3 


-0.24 


50.2 


0.081 


0.40 


1.7 


1.3 


-0.24 


59.2 


0.072 


0.45 


1.5 


1.4 


-0.26 


72.7 


0.068 




p , 

an (1450)^ TTj + 7r77 






r 

fn(1370)-» 


p 


0.265 ± O.OlSGoV 


0.270 ± 0.043GeV 


0.070 ± 0.030GeV 


0.214 ± 0.120 GeV 


0.055 ± 0.009GeV 


0.20 


0.242 


0.192 


0.119 


0.034 


0.063 


0.25 


0.250 


0.204 


0.107 


0.031 


0.057 


0.30 


0.258 


0.214 


0.104 


0.029 


0.055 


0.35 


0.273 


0.232 


0.098 


0.034 


0.058 


0.40 


0.263 


0.233 


0.098 


0.040 


0.054 


0.45 


0.272 


0.253 


0.124 


0.084 


0.056 




'"fn(1710)-.'"' + '^* 


9ao(980)KK 


9ao(980)7r77 


9/o{980)ifX 


3/o(980)7r,r 


0.137 ± 0.008GeV 


2.18 ± 0.12GeV 


1.89 ± 0.75GeV 


4.72 ± 0.82GeV 


1.12 ± 0.69GeV 


0.20 


0.177 


3.62 


2.84 


4.51 


0.31 


0.25 


0.156 


3.09 


2.34 


4.21 


0.45 


0.30 


0.140 


2.70 


1.95 


4.05 


0.61 


0.35 


0.151 


2.00 


2.06 


3.74 


0.78 


0.40 


0.129 


1.59 


2.06 


3.63 


1.00 


0.45 


0.141 


1.12 


2.13 


3.66 


1.27 



TABLE II: The results of the best fit analyses for non-derivative couphng case. The experimental data of the scalar meson 
decay widths used are cited in Ref. f^. For the 6p, we search the best fit value in the range (54.7± 18)°. The values of ga^Qm)-nr} 
etc. are ones obtained for (p ao(980)7/7r''777 and (j) ^ /o(980)7/7r''7r''7 decay analysis. 

V. CONCLUSION 

From the invariant mass distribution analysis of radiative decays (j) ~^ oq (980)7 ~^ ir^ru and (f> — > /o (980)7 ~^ tt'^tt'^j, 
we obtain the results that the KK loop diagram contribution for 0ao7 and ipfoj couplings are dominant for both 
non-derivative and derivative SPP coupling cases. We assume that cj) — > ao (980)7 ^^id (j> /o (980)7 decays are caused 
through only the KK loop diagram, and then we get the results Eqs. (18), (18'), (27), (27') of SPP coupling constants. 





non-derivative coupling 


derivative coupling 


9aKK 


2.18 ±0.12 GeV 


9.04 ±0.50 GeV-i 




1.89 ±0.75 GeV 


5.79 ± 2.32 GeV-i 


9fKK 


4.72 ± 0.82 GeV 


20.0 ±0.50 GeV-i 


9 fair T 


1.12 ±0.69 GeV 


2.43 ± 1.50 GeV^i 



We consider that the scalar ao(980) and /o(980) are qqqq states and mix with high mass scalar mesons considered as 
qq states. The low mass scalar and high mass scalar mesons are considered to mix through the inter-mixing parameter 
Aqi- In our mass formula, we obtain the mixing angle 9a, 9k and mixing parameters Rff,(g80)NN^s using the mass 
values of low mass scalar mesons (ao(980), if q (800), /o(600), /o(980)) and high mass scalar mesons and glueball 
(ao(1450), K^iim), /o(1370), /o(1500), /o(1710)). We tabulate these values for Am = (0.30 ^ 0.35) GeV^. 

e»„ = (15.0 <-> 17.8)°, = (13.8 ^ 16.4)° 
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Am (GeV^l 


A 


A' 


A" 


dp 


ao(980)— >7rT7 + ifif 




GcV^^ 


GeV^^ 




0.075 ± 0.025GeV 


0.20 


5.8 


1.2 


—0.26 


41.2 


0.093 


0.25 


8.3 


0.57 


—0.37 


36.7 


0.170 


0.30 


7.2 


0.51 


-0.33 


36.7 


0.124 


0.35 


6.0 


0.54 


-0.35 


36.7 


0.068 


0.40 


5.1 


0.57 


-0.36 


36.7 


0.054 


0.45 


4.2 


0.59 


-0.34 


50.2 


0.051 




p , 

an(1450)^ ^^ + ^3 








p 


0.265 ± O.OlSGoV 


0.270 ± 0.043GeV 


0.070 ± 0.030GeV 


0.214 ± 0.120GeV 


0.055 ± 0.009GeV 


0.20 


0.276 


0.239 


0.025 


0.009 


0.057 


0.25 


0.274 


0.262 


0.045 


0.083 


0.068 


0.30 


0.279 


0.267 


0.036 


0.099 


0.053 


0.35 


0.279 


0.266 


0.028 


0.090 


0.055 


0.40 


0.279 


0.266 


0.025 


0.082 


0.059 


0.45 


0.276 


0.246 


0.022 


0.056 


0.057 




^fn(1710)-.'"' + '^'^ 


9ao(980)ifK 


£'ao(980)7rT; 


9/0(980) if if 


3/o(980)7r,r 


0.137 ± O.OOSGeV 


9.04 ± 0.50GeV~^ 


5.79 ± 2.32GeV~^ 


20.0 ± 3.48GeV~^ 


2.43 ± 1.50GeV~^ 


0.20 


0.158 


7.80 


7.10 


8.63 


0.39 


0.25 


0.106 


11.3 


9.42 


11.6 


0.65 


0.30 


0.152 


9.65 


8.02 


9.94 


0.87 


0.35 


0.152 


7.85 


6.45 


8.23 


1.03 


0.40 


0.152 


6.45 


5.24 


6.94 


1.28 


0.45 


0.132 


5.07 


5.45 


5.67 


1.39 



TABLE III: The results of the best fit analyses for derivative coupling case. The experimental data of the scalar meson decay 
widths used are cited in Ref. 5]. For the 9p, we search the best fit value in the range (54.7 ± 18)°. The values of (7a(980)7r77 etc. 
are ones obtained for (f) — + ao(980)7/7r''777 and <f) —* /o(980)7/7r''7r'^7 decay analysis. 





Inn 


Ins 


In 




fs' 




fa 


/o(600) 


-0.98 <-> -0.97 


0.05 ^ 0.07 


0.20^ 


0.23 


-0.06^ 


-0.08 


-0.00 <-> -0.01 


/o(980) 


-0.05 ^ -0.07 


-0.93^-0.90 


0.11 ^ 


0.14 


0.34 ^ 


0.40 


- -0.02 


/o(1370) 


0.13 ^ 0.16 


-0.25 ^ -0.29 


0.48 ^ 


0.49 


-0.83 ^ 


-0.80 


- 0.02 


/o(1500) 


-0.02 ^ -0.03 


-0.03 ^ -0.05 


-0.09^ 


-0.10 


-0.02 ^ 


-0.03 


- 0.99 


/o(1710) 


-0.16 ^ -0.19 


-0.25 <^ -0.30 


-0.85^ 


-0.82 


-0.44^ 


-0.43 


-0.10 ^ -0.12 



The fact that /o(980) state considered as the Jms state mainly has the rather large fs' component 
with sign opposite to the Jns one suggests a possibility that 9 f^KK / 9aoKK '^^'^ be about 2, because 

9foKK/9aoKK = l("^-R/oWS + ^'-R/oW + V^A' Rf,,s' + A" Rf^a) I {A COS 9 a - A'sin6'a)| and Oa > 0. In our 
model, /o(1500) meson is considered as glueball, and /o(1370) meson is almost fs' state with rather large /n' 
component and /o(1710) meson is almost /n' state with rather large fs' component. 

Because 9fgKK^ ^'^'^ related to the mixing parameters Rf„Ns's and coupling strengths A, A' , A" and rj-rj' mixing 
angle 9p, we executed the best fit analyses of the various SPP decays in the wide range of the A12 value for non- 
derivative and derivative coupling cases. The best fit values of A's and 9agKK'^ tabulated in Table II and III. These 
results suggest that the non-derivative coupling seems to be reasonable than the derivative one and the intcr-mixing 
parameter A12 is rather large 0.30 0.35. 



[1] N. N. Achasov and V. N. Ivanchenko, Nucl. Phys. B315, 465(1989). 

F. E. Close and N. Isgur and S. Kumano, Nucl. Phys. B389, 513(1993). 

N. N. Achasov, V. V. Gubin, and V. I. Shevchenko, Phys. Rev. D56, 203(1997). 
[2] N. N. Achasov and V. V. Gubin, Phys. Rev. D56, 4084(1997). 

R. Delbourgo, Dongsheng Liu, M. D. Scadron, Phys. Lett. B446, 332(1999). 

N. N. Achasov Nucl. Phys. A675, 279c(2000). 



16 



[3] D. Black, M. Harada and J. Schechter, Phys. Rev. D73, 054017(2006) a rXiv:hep-ph/0601052| . 
[4] M. N. Achasov et al. [SND Collaboration], Phys. Lett. B479, 53(2000). 

R. R. R. Akhmetshin et al. [CMD-2 Collaboration], Phys. Lett. B462, 371, i-bid. 380(1999). 
A. Alosio et al. [KLEO Collaboration], Phys. Lett. B536, 209(2002); Phys. Lett. B537, 21(2002). 
[5] W.-M. Yao et al. (Particle Data Group), ,L Phys. G 33, 1(2006). 
[6] D. Black et al, Phys. Rev. D 59, 074026(1999); 

D. Black, A. H. Fariborz and J. Schechter, Phys. Rev. D 61, 074001(2000). 
[7] T. Teshima, I. Kitamura and N. Morisita, J. Phys. G: Nucl. Part. Phys. 28, 1391(2002). 
[8] T. Teshima, I. Kitamura and N. Morisita, J. Phys. G 30, 663(2004). 
[9] T. Teshima, I. Kitamura and N. Morisita, Nucl. Phys. A759, 131 (2005). 
[10] R. J. Jaffe, Phys. Rev. D15, 267(1977). 

M. Alford and R. L. Jaffe, Nucl. Phys. B578, 367(2000). 
[11] M. Harada, F. Sannio and J. Schechter, Phys. Rev. D54, 1991(1996). 
[12] D. Aston et al., Nucl. Phys. B296, 493(1988). 

[13] D. V. Bugg, A. V. Sarantsev and B. S. Zou, Nucl. Phys. B471, 59(1996). 
[14] C. Amsler et al.. Euro. Phys. J. C23, 29(2002); Phys. Lett. B353, 571(1995). 

D. Barberis et al., Phys. Lett. B471, 429(2000); Phys. Lett. B479, 59(2000). 

D. V. Bugg, A. V. Sarantsev and B. S. Zou, Nucl. Phys. B471, 59(1996). 



